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^jr^. The low-energy properties of a homogeneous one-dimensional electron system are completely 

' specified by two Tomonaga-Luttinger parameters Kp and v^. In this paper we discuss microscopic 

estimates of the values of these parameters in semiconductor quantum wires that exploit their rela- 
tionship to thermodynamic properties of the quantum wire. Motivated by the recognized similarity 
OO ' between correlations in the ground state of a one-dimensional electron liquid and correlations in 

a Wigner crystal, we evaluate these thermodynamic quantities in a self-consistent Hartree-Fock 
approximation. According to our calculations, the Hartree-Fock approximation ground state is a 
Wigner crystal at all electron densities and has antiferromagnetic order that gradually evolves from 
^ ■ spin-density- wave to localized in character as the density is lowered. Our results for Kp are in good 

' agreement with weak-coupling perturbative estimates Kp'^''^ at high densities, but deviate strongly 

I , at low densities, especially when the electron-electron interaction is screened at long distances. 

^ ■ Kp°'^^ ^ n^^^ vanishes at small carrier density n whereas we conjecture that Kp -^1/2 when n 0, 

implying that Kp should pass through a minimum at an intermediate density. Observation of such 
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a non-monotonic dependence on particle density would allow to measure the range of the micro- 
scopic interaction to be determined. In the spin sector we find that the spin velocity decreases with 
increasing interaction strength or decreasing n. Strong correlation effects make it difficult to obtain 
fully consistent estimates of Va from Hartree-Fock calculations. We conjecture that Va/vF oc n/Vo 
in the limit n — > where Vo is the interaction strength. 
PACS : 71.10.Pm, 71.10.-w, 71.10.Hf, 71.45.Lr 



I. INTRODUCTION AND OVERVIEW 



^ ■ It has been knownpfor some time that one-dimensional (ID) metals are different from their higher dimensional 
' Fermi-liquid cousinsElD. It is generally believed that at low energies and long wavelengths one dimensional electron 
systems can, under very general circumstances, be described as 'Tomonaga-Luttinger (TL) liquidsO, although it 
has nearly always been difficult to provide incontrovertible experimental evidence. Interest|ii TL liquids has betai 
heightened in recent years by new j|Jwsical realizations, including quantum Hall edge systemsoij, carbon nanotubesB'Ll, 
I ■ and semiconductor quantum wireaS^EI in particular. Like Fermi liquid theory, TL theory can be used to relate low- 
, temperature, low-frequency, long-wavelength properties to a small number of parameters in which the microscopic 
"Y^ ■ physics of particular systems is encoded. For example, TL thfiory predicts that for continuous-one-channel quantum 
' wires the quantized conductance is rcnormalized by the factonlj Kp at frequencies larger thanEil vp/L ( L is the wire 
^ : length and Up the_Eermi velocity). Surprisingly, low-energy orthogonality catastrophes lead to spectral functions that 
I follow powejjJiawsEjtl specified in terms of the same parameter. In many cases (up to logarithmically slowly varying 
^ prefactorgl3llil associated with the presence of marginal operators like back scattering in spin sector) non-universal 
Q power laws specified by TL theory parameters are also predicted for the behavior of correlation functions at distances 
^ much larger than thepS«ai;ial range of interactions. (The strictly infinite range Coulomb interaction case requires 
' special consider ations.tStS). Microscopic theory still has an important role at low energies however, in estimating the 
values of these parameters. This is especially important because the distinction between Fermi liquids and Luttinger 
liquids on the basis of a set of experimental data over a limited temperature or energy range is sometimes subtle, 
?H and the range of energies over which TL behavior is expected is often not accurately known. Approximate values af 
expected TL parameters can play a role in determining whether or not an experimental result reflects TL behaviom. 
In addition, as this approximate calculation shows, the problem of understanding the value of the two independent TL 
parameters of a homogeneous one-dimensional electron system is a challenging many-body problem that is interesting 
in its how right. 

Four TL parameters characterize the low energy properties of interacting spinful electrons moving in one channel. 
For the charge (v ~ p) or spin [ly = a) sector, the parameter Ki, fixes the exponents for most of the power laws and 
Vi, is the velocity of the long wave length excitations. Symmetries in the charge or in the spin sector reduce the numbet 
of independent parameters in the case of a one-dimensional electron gas system : spin rotation invariance enforcesEH 
Kc, = 1 while Galilean invariance implies thatt^ Vp = vp /Kp . The later identity does not apply, for example, in 
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lattice models since it requires continuous translational invariance; in that case Vp and Kp must be determined 
independently. This leaves Kp and v^, as the only two independent TL parameters for single-channel semiconductor 
quantum wires since they can be accurately described by a continuum envelope function approximation. 

In Fermi liquids a traditional and successful strategy has separated the phenomenological application of Fermi 
liquid theory from the microscopic evaluation of its parameters. To date most theoretical TL activity has focused on 
phenomenological applications; confident interpretation of experiments will require reliable microscopic estimates of 
the theory's parameters for the various physical systems of current interest. The evaluation of Fermi liquid parameters 
in two and three dimensional metals is one of the classic eacJy topics in many electron physics, with considerable recent 
progress coming from quantum Monte Carlo calculationsli3. Still, useful physical insight and reasonable accuracy has 
resulted from less computationally cumbersome approaches. In this paper we discuss what can and cannot be learned 
about the values of TL parameters in semiconductor quantum wires and the physics of their dependencies on system 
geometry using unrestricted Hartree-Fock estimates of sppund state energies. The Hartree-Fock approximation can 



yield very reliable estimates to the boundary exponentsEj describing tunneling into the end of quantum wires and, 
as a microscopic approach, gives information about quantities not reliably accessible in the TL formalism, including 
absolute values for the prefactors of power laws. 

For non-interacting electrons the TL parameter -fCp-i= 1 • With repulsive interactions its value should decrease and 
go to zero in the limit of very strong or long rangingE3 interactions. For microscopic interaction potential V{x — x') 
the formula 



K-^ = v/l + 4my(fc = 0)/7r2n (1) 

is commonly used in the literature. It depends on the carrier density n, the effective mass m and the fc = Fourier 
component of the interaction 

V{k) = I dxV{x)coskx . (2) 



Relation can be motivated by lowest order perturbation theory, or by the random phase approximationS, though 
it misses the Fock contribution for spinful electrons in ID. Any naive higher order perturbative contribution is 
divergeni— .only the infinite subsums that are conveniently captured using a perturbative renormalization approach 
are finitecj). Eq. (|^) completely ignores the renormalizing influence of short wave length modes in determining the 
actual values of the effective interaction. Higher order perturbative renormalization group calculations demonstrate 
how the interaction parameters are coupled and renormalize as short wave length contributions are integrated outCJ. 
One important example is back scattering, across the Fermi line, of opposite spins Fermions, the so called gi-process, 
that spoils the separate conservation of the number of left and right moving particles of a given spin and therewith 
is not included in the TL model. This interaction, which is finite in leading order perturbation theory, scales-ito zero 
during renormalization, restoring SU(2) spin rotation symmetry and the validity of the TL at low cncrgicsEj. Even 
for a model of spinless electrons, the parameters will rescale at low energies, reflecting other irrelevant operators-that 
are omitted in the TL model such as non-linear dispersion of the kinetic energy in the microscopic HamiltoniarJj. As 
a consequence Eq. (0) cannot be used to estimate Kp when interactions are strong. f— , 

How interactions influence the spin sector is even less certain. According to textbook knowledgec^ the spin velocity 
would-Jiot be altered by interaction forces that act only on spatial coordinates and thus not in the spin sector. Other 
worktj includes exchange contributions in the Bose form of the Hamiltonian in a way violating the SU(2) invariance 
property, K„ — 1 . On the other hand, changes in v^/v-p are quite crucial to various physical properties. It influences 
for instance the magnetic susceptibility, the g-factor, and_spin transport properties. TheJatter are particularly 
important foe, potential one-dimensional spintronic devicesEEI. In one-dimepsional channelsEj for example the spin 
conductanceES, and Rashba precession in the presence of spin orbit couplin^J, depend on ■ Most directly Va can 
be measured by inelastic Raman scattering in the 'depolarized' configuration with perpendicularly polarized incident 
and outgoing lightE30. 

To date relations between the microscopic electron-electron interaction and resulting IL parameters have been 
established for modela-pf primarily theoretical interest, such as the Kondo Lattice ModelEj, the Hubbard modelES, 
and the t — J-model.Eil For the latter two models the ground state energies are known exactly, either analytically in 
certain limiting cases or by solving the Bethe-Ansatz equations numerically. For these repulsiMa short range interaction 
models Kp is found to be confined to the range 1/2 < Kp < 1 and it has been arguedE^ in the limits of either 
infinite interaction strength or vanishing particle density are equivalent to non-interacting spinless Fermions with kp 
being replaced by 2kp so that Kp ^1/2 in either of these limits in order to recover the correct asymptotic decay of 
the density-density correlation function. For the t — J-model, TL parameters away from the supecj-symmetric point 
{J/t = 2) have been obtained by using ground state energies from exact diagonalization caleulationsE3. The Sutherland 
model for spinless Fermions, where V{x) — X/x^ , has proven to be a TL at low energiesES. The asymptotic decay of 
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its one-particle Greens function implies that Kp = 2/ (1 + \/l + 2A) with 1 > Kp > for repulsive interactions, A > 0. 
The compressibility of this system is proportional to Kp and satisfies Eq. (^) below. For quantum wires with long 
but finite range Coulomb interactioas the values of the TL parameters have been determined previously by extensive 
quantum Monte Carlo calculation£3. However, the limits on the number of particles and lattice points in real space 
for which these calculations can be carried out in a reasonable time places limits on the range of particle densities 
over which accurate Monte Carlo results can be obtained. In particular the low-density regime where interactions are 
strongest is difficult to reach. In this work we also exploit the thermodynamic relations between the uniform static 
compressibility k and the TL parameter in the charge sector. For quantum wires we have (cf. Refs. [ 

1 ^ d^{Eo/L) _ nvp^ _ ■KVF_ , . 

K ~ dn^ ^ 2 Kp ^ 2 K^p ' ^ ' 

where the last equality uses Galilean invariance. _ 

Density-density correlation function calculations performed within the TL model have been usedEj to show that 
long range interactions V(|a;|) ~ \x\~°' put the one-dimensional Fermion system ground state into a Wigner crystal 
when a < 1 , irrespective of the strength of V . The Coulomb case, a = 1 , is marginal and density-density correlations 
decay extremely slowly, slower than any power law at large distances. This observation suggests that the ground state 
energy Eq and therefore any T = thermodynamic property should be accurately estimated by the uarestricted 
Hartree-Fock approximation, for which the ground state is a Wigner crystal. At high densities it is knowncJ that the 
Hartree-Fock (HF) approximation reproduces the leading order perturbative RG result for Kp. Since it spontaneously 
breaks translational symmetry, the HF approximation correlation functions have infinite range. The tails of the 
correlation functions are therefore given slightly incorrectly. Much more important for the energy however, is the 
accurate estimate of the magnitude of the short-distance correlation fimction oscillations, illustrated in Fig. ^ The 
dominant periods for charge and spin densities are in agreement with TL model calculations, which are however, not 
able to estimate the amplitude of the oscillations. In the real correlation function, these oscillations are multiplied 
by an envelope function whose decay properties are inaccessible to Hartree-Fock theory but can be calculated from 
the TL theory. The accuracy of self-consistent Hartree-Fock energy estimates, particularly at low-densities, has been 
previously been established in other interacting electron systems.E3 

We compare our results for ground state energies and compressibilities k with estimates obtaiped within the har- 
monic approximation to the classical Wigner crystal, cf. (|^,|^) and within perturbation theoryrJ. The perturbative 
expression ( p^ we use (see below) for the charge sector TL parameter Kp turns out to be surprisingly accurate over 
a wide range of carrier densities including typical experimental ones. Only below densities corresponding to values 
Ts > 1.3 of the usual rg-parameter used to measure the interaction strength in metals (see below), do we find smaller 
K in the selfconsistent HF solution than given by Eq. (|l|) and (|). At even lower densities n ^ , where R is the 
long but finite range discussed below that we use for the electron-electron interaction, Eq. ( p^ ) predicts that k goes 
to zero whereas HF theory yields a finite limiting value, after k has passed through a minimum. This minimum is 
also reproduced by the harmonic approximation. In the latter approximation, however, k diverges as n . We 
shall give arguments supporting the conjecture that the true Kp stays finite as n — > and indeed that Kp —i- 1/2 
in this limit. 

We also analyze the spin sector and compare different approaches in the attempt to determine the spin velocity . 
The simplest estimate is again low order perturbation theory for the magnetic susceptibility. Other estimates can be 
obtained by starting with the assumption that the system is close to an antiferromagnetic Heisenberg spin chain at 
low densities. It turns out that the estimates for Va- that follow from different plausible schemes difler substantially. 
Furthermore, at smaller particle densities they are not in good agreement with earlier QMC estimated. Although we 
are able to conclude that is small at low-densities, our results for the spin sector do not substantially improve on 
existing Monte Carlo results. 

II. MICROSCOPIC INTERACTION 

Now we discuss the form of realistic interactions V along the x-direction of a quantum wire. At inter-particle 
separations \x — x'\ larger than the diameter d of the wire F(|x — cc'j > d) = e'^/e\x — x'\ will be of the Coulomb 
form, irrespective of the detailed shape of the transversal potential. If the material enclosing the wire is insulating 
with dielectric constant e', the Coulomb form still holds at larger distances, but with dielectric e replaced byCj e'. We 
assume here equal dielectric constants e ~ e' , the case that applies to gatedB as well as in cleaved edge overgrowth 
structureaj. Eventually, at an inter-particle separation exceeding the distance R to the closest metallic system that 
is extended in the direction along the wire. This metallic screening can be supplied by carriers in nearby metallic 
gates, including those used to define the quantum wire. Assuming that the screening plane and the quantum wire are 
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parallel V{\x — x'\ > R) l/\x — x'\'^ because of the formation of dipoles from image charges. The interaction (Q) 
below accounts for this cut-off at large particle separations. 

At distances \x\ ^ d shorter than the wire width, the precise transverse form of the electronic wave function 
influences . For example, in 2D hetero-structures all electrons share a common growth direction wave function. 

If this wave function is approximated by aJiarmonic oscillator with characteristic length d, and thickness in the 
2nd confined direction is neglected, we haveEl V^nix) = {e^l2^ed)e 1^'^ Ko{x'^ /M^) where Kq denotes a Bessel 
function. It is more realistic to include finite thickness in both confined directions. For example, for 3D win 
with circular cross sections, a model which might be appropriate for cleaved-edge-overgrowth systems, we havd 
V,u{x) = (eVV2Aed)e"'/2'^'erfc(a;/V2d) again using harmonic oscillator ground states of widths d/\/2 , now for 
both of the transverse directions. Neither of these forms is smooth at x = , an artifact of assuming factorized wave 
functions. At short distances, corresponding to high energies, the factorization assumption needs to be refined. It 
leads to the unphysically slow decay of the Fourier transforms V^-oik) = {e^ /evF)e''''^~ Koik^d"^ /2) , and Vanik) — 

, as seen in Figure ^ Ei is the exponential integral. Also included in Figure |^ is the more 

realistic form, 

V{k) = IXoikd) ~ KoikVd^ +4R^)\ , (4) 

which accounts for the image potential term from a remote screening plane separated by i? 3> d. Its real space form 



is displayed in Eq. (A4). This interaction remains finite for k and decreases more quickly for large k and will 
be used in the present work. Here, and in the following we measure the interaction in units of the Fermi velocity so 
that V is dimensionless in Eq. (^. Its strength, in comparison with the kinetic energy, scales with the dimensionless 
parameter :— 1/ (2naB) — vr/ (4fcFaB) depending on density n and the effective Bohr radius ob . In our calculations 
we assume ob = 2(i . 
Thus, two parameters 

R/d and kpd (5) 

characterize the range and the strength of our model interaction (|^), respectively. They both can be extracted quite 
reliably from experiment, R from the sample lay out and d from the energy ~ l/md^ of inter-subband excitations. 
Typical distances to metallic gates and typical wire width, as reported eg. in [ |^ and [ ^, correspond to values for 
R/d ranging from 5 to 14. Typical single wall carbon nanotube systems, on the other hand, would correspond to 
much larger R/d values because of their extremely small diameters. Many of our calculations are for R/d = 5.66 or 
35.36. Note also that electron densities should be sufficiently low ( kpd < \f2 within the parabolic approximation for 
the transverse confinement) to prevent occupation of the second subband. 



III. GROUND STATE ENERGY 

According to (^) we need to calculate ground state energies for different particle densities. In this work we employ 
the unrestricted Hartree-Fock approximation where details are described in Appendix Results of these calculations 
are included in subsequent Figures. 

The close proximity to the Wigner crystal and the Bose character of all of the low energy excitations suggests 
comparing with the ground state energy density of the harmonic crystal in ID, 

ET = Et'^'^^' + - — Lo{k) . (6) 
2 27r 

Here, 

is the classical contribution and the zero point energy follows from the phonon dispersion 

u;2(fc) = - E y"(j7r/2fcF)(l - cos(jfc7r/2A;F)) (8) 
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of harmonic excitations. The primes denote derivatives w.r.t. the argument. Both, jrjciassicai j^wc pj-Qvide 

rigorous lower bounds to the true ground state energy since the quartic term of the Coulomb interaction is positive 
when expanded in a power series and the Fermionic antisymmetry constraint, ignored by Eq. (^), increases the true 
Fermionic energy further. This latter observation remains true also for spin carrying electrons since spin cannot 
provide complete antisymmetry for symmetric spatial wave functions for more than two particles. 
Figures || also include the lowest order perturbation theory estimate 



obtained by taking the Hamiltonian's expectation value in the non-interacting electron state, to obtain the positive 
Hartree, and the negative exchange contribution. The variational principle ensures that Eq. is a rigorous upper 
bound to the ground state energy. The true ground state energy must lie between these two bounds. 
The energy densities are plotted in dimensionless units: 

which has the value I/Stt without interactions. The HF-energies e[^^, seen in Figure || a for Rjd = 5.66 and in 
Figure |^ b for = 35.36 , approach this value in the weakly interacting high density limit, fcp — > 00 . 

For densities above fcpd ^ 0.5 , corresponding to ^ 0.8 , Bq ^ agrees quite well with the perturbative estimate 
. This is despite the fact that the self consistent charge density modulation already shows significant amplitude in 
this regime as seen in Figure ^ below. 

Below kjfd ^ 0.3 (rg ^ 1.3), the three approximations start to spread apart significantly. As /cp perturbation 

theory result diverges, Cg*"^* '-^ 1/fcp , while e^'^ ~ kp^ goes to zero. The mean field result e^^ approaches a finite 
value, a result which seems plausible since the zero point quantum fluctuation energy exceeds the classical interaction 
energy of the Wigner crystal at low densities for interactions decaying faster than ^ l/\x — x'\'^ at large particle 
separations. We speculate that realistic quantum wires, which never have strictly infinite range interactions, always 
cross over into the hard sphere gases at sufficiently low densities. For this latter system it is known that the ground 
state energy approaches cq 4/37r « 0.4244 (cf. [ ||J4|]), irrespective of the particle type, Fermionic or Bosonic, 
and irrespective of the particle's spin. The 'radius' of the hard spheres is unimportant when hp —^ . Among the 
approximations discussed above e^^ is the only one that stays finite in this limit, though the limit it approaches is 
larger than 4/37r . 



IV. TL PARAMETER Kp 

Figure ^ shows 

1 /2 

l/Kp = (|[fc2e;;(fcF) + 6(A:pe^(fcF) + eo(fcp))]) (11) 

versus kwd for R/d — 5.66 (Figure^ a) and for R/d = 35.36 (Figure ^jb). Eq. ( [Tl| ) follows from relation ^ together 
with (y|); the primes again denote derivatives w.r.t. the arguments. Also included in these figures is the result from 
expression (|l|) and the perturbative estimate 

l/XP-t = (l + (21/(fc^0)-F(fc = 2fcp))/7r)i/2 , (12) 

which follows from Eqs. (^) and (||)). Note that only this form, with the Fock term included, satisfies the physical 
requirement that spinless Fermions cannot feel contact interactions and that therefore Kp equals unity for this model. 
Since the factor of two in ( |l^ ) is absent in the spinless case, Eq. ([l^ ) indeed fulfills this Pauli principle requirement, 
unlike Eq. (|l|) . Figure ^ also includes 1 /K"// and 1 /KJ"^ , calculated from the corresponding harmonic crystal energy 
estimates of Eqs. (0) and (||) using Eq. (||). 

Over a wide range of densities, including the typical experimental regime, all of these approximations coarsely 
agree, though none of them can provide a rigorous bound on the exact compressibility or Kp . As for the ground 
state energies, the approximations start to deviate severely from one another at smaller densities, corresponding to 
fg ^ 1.5 . Both HF and harmonic estimates-show non-monotonic behavior of ^/Kp as a function of density, in 
agreement with recent quantum Monte CarloEII calculations. If, as we have conjectured, cq approaches a constant 
for fcp — > , 



5 



i^p(fcF ^ 0) = (3^eo(fcF -> 0))-i/2 . (13) 

Note that the HF compressibihty approaches a constant in the low density limit. Conjecturing again that at mean 
particle separations exceeding the interaction range, kp <C n/2R , the system crosses over into the hard core Bose gas 
with eo 4/37r Eq. (|l|) would yield 

K,~^l/2 . (14) 

We note, for example, that the Hubbard model approaches ( [T^ ) at snaalL fillings, independent of the interaction 
strength U. The same holds true for the Fermi gas with contact repulsionl3'C3. It is an important observation that the 
limiting value K^^{kp ^ 0) w 0^29 to 0.35 for R/d = 50 . . .8 clearly exceeds 1/8 which would be the limiting value 
for the extended Hubbard modelc3 that has both on-site and near-neighbor interactions. On the other hand, as seen 
in Figure ^, the minimum value for Kp at about kpR ~ 1 is considerably smaller than 1/2, so that, contrary to the 
Hubbard model, the limit (|lj) would have to be approached from below with decreasing carrier densities in quantum 
wires. 

Upon inspecting Figures |^ more closely a regime can be identified at densities somewhat above the maximum of 
1 /K^^ where 1 /K^^ exceeds 1 / ifP'^''* . As seen in Figure ||, the relative increase in stiffness appears along with the 
occurrence of significant 4fcF-periodic contribution to the charge density modulation. Figure |^ shows the two lowest 
Fourier coefficients 

Qj ^ Q^{q^ j'^kp) = {-lygiiq = j2fcF) = Q-j (15) 

for j = 1,2 in units of the mean density go — 2kpl'K . In view of (p^), which follows from Eq. (A7) of the Appendix, 
4A;F-periodic modulations of the charge density (x) + qi {x) are given by the 3=2 contribution in Figure ||. The 
appearance of a substantial j = 2 Fourier component, at k-pd ~ 0.5 , marks the crossover from spin-density-wave to 
Wigner crystal self-consistent solutions of the HF equations. A similar conclusion has been drawn from the, extremely 
slow spatial decay of the density-density correlation function in the presence of long range interactionsEj and from 
recent quantum Monte Carlo studied^. With smaller 1/i? this regime of Wigner crystal- like states marked by 
enhanced stiffness extends down to smaller densities and becomes more pronounced. The variation of XjK^^ with 
R is depicted in Figure ^ for the density kyd = 0.15 . At /cFi? S> 1 all of the approximate estimates are consistent 
with the logarithmic increase l/ifp ~ -^/log Rjd suggested by perturbation theory. For kpR^ 1 , the electrostatic 
energy is so dominant that the energy and Kp are relatively insensitive to correlations. 



V. SPIN SECTOR 



As mentioned already in the introduction, it is much more difficult to estimate how interactions influence the spin 
sector, and particularly its low energy TLj-parameter Va, than it is to estimate the charge sector parameter. In 
the model originally proposed by Luttingera with left and right going particles treated as distinguishable, the spin 
velocity is unrenormalizedKI, Va = vp, because the exchange term vanishes leaving magnetic properties-pf the system 
independent of interactions. For the Hubbard model, on the other hand, is known that the spin velocitycZl is dependent 
on interactions and particle density, n , vanishing like at small density n for any finiteJuateraction strength. The 
spin TL parameter is related to a thermodynamic quantity, the magnetic susceptibility, byCj 

2K 

X^4(^V9^eo(m))-i = ^ (16) 

where m = (n^ — ni)/n is the magnetization per particle and eo is the dimensionless ground state energy density 
as defiiied in Eq. dlQ). Relation (llfl) actually holds for any interacting electron system in the single channel TL 
phasdfj'LJ. Evaluating eQ{m) for the microscopic model perturbatively would give 

ijvp = 1 - V{2kp)/n , (17) 

using (^|) and K„ = 1 . Alternatively one also could impose a spin current i — (riai — n^i — ^^lT + ^^lI ) / n per particle 
[n-R/Ls are right/left moving densities of spin s) and measure the change in ground state energy 

X, = 4{7r'vFd^,eo{i))-' = . (18) 

nVrrK„ 
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Perturbatively this gives an unchanged spin velocity, a result that simply reflects the fact that lowest order perturbation 
theory cannot describe drag effectsEj between the density fluctuations of opposite spins and thus leaves the system 
Galilean invariant in spin sector. Solving equations ( |l6| ) and (|l^ for K^r yields K^r > 1 as a perturbative result for 
repulsive interactions which would contradict SU(2) invariance in a TL model. To enforce the SU(2) symmetry we 
can combine equations (Iq) and (Uq) and solve for by eliminating K„. The result, 



= -{xXiT^'^ - vyJl-V{2k^)/n, (19) 

TT ' 

indeed agrees clearly better with the QMC data0 than (|l^. Eq. (|l|) is included in Figure |. At V{2k-p) = vr the 
perturbative estimate uP°''* vanishes and for smaller fcp the Fock term in ( [T9| ) favors a spin polarized ground state. 
This result contradicts very general arguments that guarantee a non-magnetic ground state for any non-singular 
pair interaction potential in one dimensioned. The true spin velocity should stay positive and approach zero only at 
vanishing particle density. 

That the extraction of spin velocities from HF calculations is less reliable than the extraction of charge TL pa- 
rameters, is already clear because of the incorrectly broken spin-rotational invariance in the HF ground state In the 
HF spin-density wave state we evaluate the spin-susceptibility by polarizing spins along the quantization axes. We 
can consider only cases with rational ratios of the spin-up and spin-down carrier densities. Because the periods of 
the spin-up and spin-down density waves differ in these solutions, it is more convenient to use a real space basis, 
discretizing space ip{x = Xi) —^ Tp(i) as described in Appendix^. Self consistent solutions are shown in Figure]^. 
At finite magnetization this-Sjtructure contains now 'defects', reflecting the loss of the 4fcF-periodic component in the 
charge density modulations^. At least N — electrons have been considered on M = 401 grid points, the smaller 
particle densities are based on iV = 84 and M = 801 to avoid lattice artifacts to a high accuracy. These sizes 
are clearly beyond what presently can be treated with numerical many body approachesjike quantum Monte Carlo, 
but pose no problem here. Spin velocities obtained from E^^{m)/L by virtue of Eq. ( p^ ) are included in Figure ||. 
Below kpd = 0.2 it is very difficult to extract positive spin velocities. We see that selfconsistency pushes the point 
of vanishing spin velocity and the (erroneous) transition into a ferromagnetic ground state down to smaller densities 
compared to the perturbative estimate in Eq. (p^, but the transition still occurs. 

An alternative attempt to estimate the spin velocity starts from the argument that the electron spin sector would 
evolve at low particle densities towards that of an antiferromagnetic Heisenberg spin chain, as suggested by the 
staggered spin density profile found in the mean field solution. Fig |^. This argument is also suggested by the 
pronounced antiferromagnetic correlations fouad in the TL liquid spin sector, particularly for long range interactionscJ, 
and in finite pieces of one-dimensional wiresEj. Antiferromagnetic spin chains are known to represent microscopic 
models, such as the Hubbard model, at low energies and have been. iutcpsively investigated for instance by employing 
manifestly SU(2) spin rotation invariant non-Abelian bosonizationliJO. 

In the antiferromagnetic Heisenberg chain spin excitations (magnons) move at velocity 

Va = 7r^J/4fcF 

where J is the nearest neighbor coupling constant. One possibility to guess the magnitude of J is to compare the 
HF estimates for the ground state energy e^^ of unpolarized electrons with the ground state energy Cq"' of foJly 
spin polarized electrons. For the antiferromagnetic Heisenberg chain this energy difference, J(l + In 2) per spin,L3 is 
known exactly. Equating the energy differences gives 

T^pfcp HF) (20) 

1 + In 2 2 ^ ° ° ^ ^ ' 

from which 

«>F = 7r3(eP°'-er)/8(l + ln2) (21) 

follows. Eq. ( |2l| ) is included in Figure ^. The transition into the spin polarized ground state occurs at kpd — 0.19 
{vg — 2.07) for R/d = 5.66. Equation (EO) can be checked for consistency in the non-interacting limit, kp ^ oo, where 
v„ ^ vp . Magnons would move at velocity vp if J = Avpkp/n"^ = OAlvpkp . On the other hand, (eg° — e^^) — > I/tt 



in this limit so that (20) yields J — 0.30 uf^f ■ In view of the fact that the weak interaction limit is poorly described 



by the antiferromagnetic spin chain this picture seems amazingly consistent. 
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Recently, Calmels and Gold have calculated magnetic susceptibilities of quantum wires 
microscopic interaction, using standard heuristic approximations from electron gas theoryEll. By virtue of Eq. ( p^ ) 
these data allow us to extract spin velocities which turn out to be slightly larger than our v^^ values. Note 



though for a different 
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that the perturbative estimate shown in Figure 2 of Ref. [ Q uses Eq. (|T^) while in Figure || Eq. (|9|) is included. 
Compared with the data are smaller than ?;P°''* at Is^ge densities, like the Hartree-Fock data. This latter 

property does not agree with the behavior obtained using QMCEZi where v^^^ instead exceeds t'P'"'' . We conclude 
that HF and other approximations of the mean field type can provide only a qualitative guideline to . All of these 
attempts, however, agree in predicting spin velocities that depend on the interaction and decrease with increasing 
interaction strength. This result calls attention to thc-Jrequcnt assumptions in the literature that interactions not 
explicitly depending on spin would leave Va unchangedcJ. This result should show up in current experiments, such as 
those described in Ref. [|9|, where typical values for kpd w 0.3 are in the regime investigated here. As already pointed 
out in the introduction, this parameter should influence-measurable quantities, sucb_as the spin-splitting enhancement 
factor, Raman scattering in depolarized configurationEj, spin transport propertiesO, and Rashba precessionE3. 

Let us now discuss the low density limit using our conjecture that quantum wires become equivalent to the on-site 
Hubbard model in the limit of small particle densities. Hubbard model (lattice constant a) parameters, t = v-p/2kpa'^ 
and U — V{k — 0)vF/a, can be related to microscopic parameters by equating the effective mass and the Coulomb 
barrier for a tvpielectron exchange, V{k) is defined in Eq.(||). To leading order in t/U the spin velocity of the Hubbard 
model v^^ isEl 



( sin4fcFa 



so that 



U \ ikpa 



ATT/3V{k = 0) (22) 



Vp 

for small kpa , where the lattice constant is irrelevant. Note that V ^ v^^ and thus cx kp . This result agrees 
with the strong interaction limit of the continuum version, the electron gas with repulsive contact interactionsEj. 
Restoring quantum wire parameters Eq. (|2^) translates into 

fcp-o 27r fepQE 

i^F ^ 3 ln(2i?/d) ^ ' 



for R/d ^ 1 . The available QMC data are consistent with (23), though, as in the charge sector, they are not 
conclusive enough to really confirm the low density equivalence. 

VI. SUMMARY AND DISCUSSION 

Many non-trivial theoretical predictions exist in the literature on various measurable low energy properties of single 
channel quantum wires that are based on the TL model which describes the non- Fermi liquid behavior of these systems. 
Much in the spirit of the Landau theory of Fermi liquids all these predictions depend only on few phenomenological 
parameters. In the case of quantum wires with only one subband occupied there is only one parameter per degree 
of freedom, Kp for the charge sector and v„ for the spin sector. In the absence of interactions these parameters 
assume the values Kp = 1 and = vp . In this work we have investigated how the microscopic pair potential 
V{x — x') changes Kp and Wo- • We have considered a realistic, tractable, and sufficiently general form for V{x — x') , 



Eq. ( A4), that depends on the diameter d of the quantum wire which is measurable through the subband energy, and 
the distance to the nearest metallic gates R as given by the sample lay out. Our approach is to relate the two TL 
parameters to thermodynamic quantities which we estimate on the basis of self consistent, unrestricted Hartree-Fock 
(HF) approximations for the ground state energy. 

In the charge sector this strategy is found to yield reasonably accurate results. At densities corresponding to 7's 1.3 
we confirm applicability of the pert™ iative formula ( p^ for Kp . This regime includes most of the experiments based 
on semiconducting heterostructuresElEj. At somewhat smaller densities Eq. (|l2|) even overestimates Kp . In this 
regime we find enhanced stiffness compared to perturbation theory, signaling the close proximity of a Wigner-crystal 
state. For this reason quantitative corrections to Eq. ( [T^ ) may arise when R d , for example in carbon nanotubes. 
In quantum wires fabricated on the basis of semiconducting hetero-structures with gates, the perturbative formula 
may be used even down to densities 2kY/TT fa 1/R . The proximity to the Wigner-crystal state competes with the 
finite interaction range in these systems. Irrelevant or marginal operators in the microscopic Hamiltonian, such as 
non-linear single particle dispersion or backward scattering in the spin sector turn out to be unexpectedly inefficient 
to renormalize the TL parameters in the charge sector up to moderate interaction strengths. With decreasing density 
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the values for Kp clearly fall short of 1 /8 which is the minimum assumed by the extended Hubbard model including 
repulsions on neighboring lattice sites and often is considered to emulate models of finite interaction range. 

At smaller densities, however, does i^TP"''* become poor. In particular, it does not reproduce the non-monotonic 
behavior of Kp found in the HF approximation with a minimum as a function of density. Eventually, as /cp — *■ 
we conjecture that quantum wires approach the universality class of the Hubbard model with only on-site repulsion 
and that Kp -^1/2 in that limit, though, unlike the Hubbard model, this limiting value should be approached from 
below as the particle density is lowered. i-. 

The non-monotonic dependence of Kp on density predicted here should show up in any of the power lawaj revealed 
by pseudo gaps in the density of states v. Examples include the current for tunneling into the end ^ cj^^/^p"^^/^) 

or into the middle {v{uj) ~ tt;(-^p+i/^p-2)/4-j ^ single mode wire (assuming K^ — 1), and the current liy) ~ y^/^p 
flowing through a single tunnel barrier along the wire at small voltages V . Experimental observation of this non- 
monotonic dependence of the exponent would give direct experimental access to the microscopic range of the electron- 
electron interaction; the position and the height of the maximum in K~^ both depend on R . 

Our approach is less successful in estimating the spin sector TL parameter , at least when it differs considerably 
from Dp. We have discussed perturbation theory and tried to obtain meaningful estimates for Vu from the HF 
spin density wave states. The resemblance to the antiferromagnetic Heisenberg spin chain, suggested by correlation 
function considerations, suggests that exchange coupling strengths, and therefore spin velocities, can be estimated by 
comparing the ground state energies of unpolarized and fully spin polarized electrons. None of these variants lead 
to results of the same quantitative reliability as those obtained from K^^ . Conjecturing again a cross over into the 
universality class of the Hubbard model in the limit of /cp ^ yields the prediction of a linear dependence of the 
relative spin velocity Va-/v-p cx k-p/Vo on the particle density. Vq is the zeroth Fourier component of the interaction 
V{x-x') . 

It is important to knew the spin velocities for attempts to realize 'spintronic' .devices where spins rather than 
charges are txansportedcj, using for example the Rashba spin, precession mechanismEa through quasi one-dimensional 
constrictionfO. Here, in agreement with QMC estimatesEZi, we have collected strong evidence that spin density 
excitations move at speeds considerably slower than the Fermi velocity, already in present day deviceaa, where kpd « 
0.3. 
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APPENDIX A: HARTREE-FOCK THEORY 



To formulate the mean field theory we introduce single particle wave functions ip solving the Schrodinger equation 
|-^52-|-^l/H(x)|^fc,(a;)-^ dx'V,''ix,x')4,ksix')^eks^ks{^) (Al) 



with the Hartree 

2n 



= / da;' V{x-x') / dfc |^fe.-(x')P (A2) 



and the non-local exchange 



Vf{x,x') = ^V{x x') r dfc rks{^')^PM (A3) 



27r 



potentials, which itself depend on ip and thus have to be obtained selfconsistently; s , | are spin quantum numbers. 
Occupation of only the lowest subband is assumed together with periodic boundary conditions for the wire of length 
L . Parabolic dispersion for the kinetic energy in ( |Al[ ) is described by a band mass m — kp /v-p . The kinetic energy 
is not linearized. 
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For some of our calculations, particularly those focussing on properties of the spin sector, cf. Section ^ we solved 
the HF equations (Al) directly in real space, using 



V{\x\)^-{ 



1 



(A4) 



in (A^,A3) and a lattice grid of at least 401 points. Any of the results for the charge sector can be obtained either 
using a real space basis as also, slightly more efficiently, in A; -space, introduced now. Expanding 



(A5) 



into Bloch waves, and similarly the periodic potentials (A2) and (A3), yields HF equations for the coefficients Uj,k,s 



= 



2 Kp kpvp 



Uj.k,s + 



2kpTT 



3 3 



k,s I dk' 



(A6) 



^-j+j'+j" -k' ,s '^j'.k'.s 



For each s = ±1 and fc = — fcp . . . fcp inside the Brillouin zone this is an eigenvalue equation for matrices indexed by 
the band indices j which, however, inside the curly bracket, depends on the solution of (A6). Within the 'unrestricted' 
HF scheme we allow for charge and spin density wave solutions breaking the symmetry of continuous translations and 
thereby lower the ground state energy. Solutions are found to show 4fcF-periodic oscillations of the charge density 
g{x) = Qiix) + Qiix) where 



27r 

dfc \ipk^{x)\^ = + 



(A7) 



We solved Eq. (A6) iteratively, starting with a sinusoidal spin density wave u 

i M ^ '^J.o/\/2 + S(5|j|,i/2 . The final 



solution always obeys Mj,fc,T = (~l)"'%,fc.j. which in view of (A?) yields 2fcF-periodic modulations of the spin density 
{x) — Qi {x) . A typical density modulation at str ong er interaction is shown in Figure ^. 
The single particle energies £ks , obtained with (A6), determine the ground state energy 



L 



E 



dfc 

2^ 



£ks 



-y 

2 ^ 



dfc 



dxV^{x)\^ks{x)\'^ 



(A8) 



E 



fcp 



dfc 

2^ 



dx 



dx'vf{x,x')rks{x),kUx') 



Half of the interaction has to be subtracted to repair for its double counting in (A6). Differentiating Eq/L twice 
w.r.t. fcp yields our estimate for Kp , according to (||). 

Most of the results are obtained for 82 fc-points in the Brillouin zone (in some cases of very small densities we 
increased this number to 234). The Milne-rule, being accurate to 7-th order in the spacing between fc-points, is used 
for the fc-integrations. We included j = —3, . . . , 3 bands though in most cases j = —2, . . . , 2 would have sufficed 
due to the rapid decay of the Coulomb interaction in fc-space. 



APPENDIX B: HARTREE-FOCK FOR SUSCEPTIBILITIES 



In spin space we us e a lattice representation ip{x — Xi) — > '^(i) of the Hamiltonian. The M x M matrices 
Hiji representing Eq. (Al), contain contributions from the kinetic energy Ha = 2{M / TTN^vpkp and Hii±i = 

— {M /irNYvpkp , the (local) Hartree-term Ha — J2kjs ^il^^ j\E/M)\ipks{j)\'^ , and the (non-local) Fock-term H^j = 

— J2k ^(1* ~ j\L / M)tpks{i)^ks{j) , the latter acting only on spin-s wave functions. Here, V{x) = J^fL-oo + 
accounts for periodic boundary conditions {V{x) is defined in (A4)) and the real and normalized eigenvectors i^ks{j) 
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of Hij are indexed by their spin s = ±1 and momentum — fcps < k < kps = (1 + srnjk-p with k being an integer 
multiple of 2it/L and m the m agn etization per particle. The Hartree-Fock approximation to the ground state energy 
Eq then is obtained as in Eq. (Rq). 
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Figures 



FIG. 1. 

Figure 1: Charge densities n|(a;) (solid) and ni{x) (dashed) as a function of position along the wire x in units of the mean 
electron spacing for kpd — 0.15 and R/d = 5.66 . We argue that these charge densities in the broken symmetry Hartree-Fock 
states are a good approximation to typical configurations in the fluctuating one-dimensional electron liquid that does not have 
broken translational symmetry. 

FIG. 2. 

Figure 2: Fourier transforms of three difi'erent forms of microscopic electron-electron interactions, V{k) — V{k)/vF, as described 
in the text. Solid: the form (^ we use in the present work, dashed: the 2D heterostructures model, dash dotted: the 3D 
cylindrical case intended for cleaved edge overgrowth systems 



FIG. 3. 

Figure 3: Ground state energy densities eo(fcF) — {Eo/ L) /{kp/m) for R/d — 5.66 (a) and for R/d = 35.36 (b) in HF 
(solid). The perturbative estimate (Eq. (^, long-dashed) establishes an upper bound while the harmonic chain estimates, both, 
omitting (Eq. (^, dotted) or including quantum fluctuations (Eq. (^), dashed) both establish lower bounds to the true ground 
state energy. 



FIG. 4. 

Figure 4: l/Kp versus kpd for R/d — 5.66 (a) and for R/d = 35.36 (b). The same approximations are included as in 
Figure ^ together with the commonly used formula (|l|), dash-dotted. 



FIG. 5. 

Figure 5: 1/Kp versus R/d for kpd — 0.15 in HF (solid), perturbation theory (Eq. (|l^, long-dashed), and for the harmonic 
chain (Eq. (Q), dashed) and (Eq. (||), dotted). 

FIG. 6. 

Figure 6: Amplitudes for the 2fcp and Ifcp -periodic components of the charge density modulations pi and p2 versus kpd for 
R/d = 35.36 in units of po = 2fcp/7r. 

FIG. 7. 

Figure 7: Charge densities ni{x) (solid) and ni{x) (dashed) along the wire x in units of the mean electron spacing, as in 
Figure |l]but for flnite magnetization m = 1/6 per spin. 

FIG. 8. 

Figure 8: Estimates to spin velocities v^/vy , based on the selfconsistent Hartree-Fock solution : 'HF', on perturbation theory, 
Eq. (|l9[) : 'pert', and on the comparison with the antiferromagnetic Heisenberg model, Eq. (Ell) : 'J'. 
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